Simulating typical entanglement with many-body Hamiltonian dynamics 
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We study the time evolution of the amount of entanglement generated by one dimensional spin- 
1/2 Ising-type Hamiltonians composed of many-body interactions. We investigate sets of states 
randomly selected during the time evolution generated by several types of time-independent Hamil- 
tonians by analyzing the distributions of the amount of entanglement of the sets. We compare 
such entanglement distributions with that of typical entanglement, entanglement of a set of states 
randomly selected from a Hilbert space with respect to the unitarily invariant measure. We show 
that the entanglement distribution obtained by a time-independent Hamiltonian can simulate the 
average and standard deviation of the typical entanglement, if the Hamiltonian contains suitable 
many-body interactions. We also show that the time required to achieve such a distribution is 
polynomial in the system size for certain types of Hamiltonians. 

PACS numbers: 03.67.Mn, 03.67.Bg, 75.10.Pq 



I. INTRODUCTION 

Random states are a set of pure states uniformly dis- 
tributed in a Hilbert space with respect to the unitar- 
ily invariant measure. In quantum information, random 
states are used in several quantum protocols, for example, 
remote preparation of quantum states [l[ , quantum data 
hiding [2|, |3[ and quantum one-time pads In physics, 
the properties of random states arc studied in relation to 
a derivation of a canonical distribution in quantum sta- 
tistical mechanics 043 and information leakage by the 
evaporation of black holes @. One way to character- 
ize a set of states is by analyzing the distribution of the 
amount of entanglement of the states, often called the 
entanglement distribution. Recently, it has been shown 
that the entanglement distribution of random states con- 
centrates just below the maximum in several measures of 
entanglement |{|. Such characteristic appearing in ran- 
dom states is referred to as typical entanglement. 

Generation of exact random states requires exponen- 
tial resources, in the sense that a quantum circuit that 
outputs such states for arbitrary input states must con- 
tain an exponential number of elementary gates. But 
exact randomness is not always necessary for all appli- 
cations. This leads to the study of i-designs, where sets 
of states that can effectively simulate random states up 
to their i-th order statistical moments are investigated. 
Realizin g su ch sets in pol yno mial time has been studied 
in Refs. fl0l - [l5| . In Rcfs. p^ - [l4| . the average amount of 
entanglement is used as an indicator of the realizability 
of a 2-design, and it is shown that 2-designs are implc- 
mentable in polynomial time by random quantum cir- 
cuits [l5| ]., or by a simple measurement procedure on a 
weighted graph state fl4j |. 

However, when we analyze random states for under- 
standing the foundation of several models in physics 0- 
[g], random quantum circuits and weighted graph states 
are rather artificial, since their implementations in phys- 
ical systems require fine controls of many parameters in 



time. It is natural to ask if random states are approxi- 
mately realizable by a set of states randomly selected dur- 
ing time evolution governed by fundamental equations of 
motion such as the Schrodinger equation, or a master 
equation. This question has been addressed mainly in 
the field of quantum chaos by investigating the time evo- 
lution of the amount of entanglement that is referred to 
as entanglement evolution. It has been shown that the 
chaotic dynamics described by a time-dependent Hamil- 
tonian can generate the same amount of entanglement as 
that of the typical entanglement on average jig - tig . 

Is time-dependent control of parameters necessary for 
simulating the typical entanglement of random states? 
In this paper, to investigate this question, we consider 
a set of states randomly selected during the time evo- 
lution generated by time-independent Hamiltonians and 
investigate the properties of the Hamiltonians for simu- 
lating the typical entanglement. It is known that Hamil- 
tonians composed of local interactions acting on a fixed 
number of consecutive particles cannot generate the same 
amount of entanglement as that of the typical entangle- 
ment in finite time from separable initial states 0, [l!| . 
For this reason, we have to consider many-body interac- 
tions. In particular, we deal with one-dimensional spin- 
1/2 Ising-type Hamiltonians including various kinds of 
many-body interactions. We choose Ising-type Hamilto- 
nians, since the simplicity of Ising-type models allows us 
to analytically investigate the enhancement of entangle- 
ment generation due to the many-body interactions. We 
note that recently, realizations of many-body interaction 
Hamiltonians have been studied, for example, using Ry- 
dberg atoms [2fj| . 

Since the time-independent Hamiltonian dynamics 
change only the phases of states in the eigenbasis, it 
cannot generate exact random states. Our interest is 
how well the time evolution generated by the time- 
independent many-body interaction Hamiltonians ap- 
proximates the properties of the typical entanglement of 
random states. Due to the extremely high concentration 
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of the typical entanglement, it is necessary to investigate 
Hamiltonians which can generate high entanglement on 
average. The Ising-type Hamiltonians have the possibil- 
ity to generate even larger amount of entanglement than 
that of the typical entanglement [23[. This is another 
reason why we consider the Ising-type many-body inter- 
action Hamiltonian. 

In this paper, by examining the entanglement evolu- 
tion for several kinds of many-body Ising-type interac- 
tions, we show that many-body interactions can dramat- 
ically enhance entanglement generation. We also show 
that the average and standard deviation of the entan- 
glement distribution can be comparable to those of the 
typical entanglement in polynomial time for certain types 
of many-body interactions. 

This paper is organized as follows. In Section [Hi 
the measure of entanglement and the Hamiltonians in- 
vestigated in this paper are presented. In Section IIII1 
we present the results of the entanglement evolution for 
Hamiltonians composed of neighboring n-body interac- 
tions. In Section IIV1 we present the results of the en- 
tanglement evolution for Hamiltonians composed up to 
neighboring n-body interactions. We compare the en- 
tanglement distributions of our models and the typical 
entanglement of random states in Section [Vj We present 
a summary in Section [VII 

II. PRELIMINARY 

A. Measure of entanglement 

The linear entropy of reduced density matrices is often 
used for analyzing the rcalizability of approximate ran- 
dom states, since it is a good indicator of the 2-design 
[l2T - [l4j ]. Accordingly, we use the average linear entropy 
of one-site reduced density matrices, known as the Meyer- 
Wallach measure of entanglement. For a pure state |<I>) of 
a iV-spin system, the Mayer- Wallach measure is denned 
by 

2 N 

where pi Tr_,j |$) ($| is the reduced density matrix 
at the i-th spin, and Sl(p) := 1 — Trp 2 is a linear en- 
tropy [U [22j ■ (The partial trace Tr^i is taken for the 
degrees of freedom all spins except for the i-th spin.) 

The Mayer- Wallach measure satisfies < Emw(\®)) < 
1. Emw(\$)) = if and only if |$) is a separable state, 
and Emw(\&)) = 1 if and only if |$) is local unitarily 
equivalent to a GHZ state. For the typical entanglement 
of random states, the entanglement distribution accord- 
ing to the Meyer- Wallach measure is obtained in Ref. 
[lq . Its average and standard deviation are given by 

3 

( £ Mw)mnd = 1 _ 2 N + 1' ^ 



and 

CTrand = O ( ^ ) . (2) 



B. Generalized Ising model 

We present generalizations of a one-dimensional spin- 
1/2 Ising Hamiltonian for many-body interactions. We 
consider a one dimensional A-spin system, in which each 
spin is specified by the index of its site i, and interac- 
tions of spins act only on consecutive spins. We impose 
a periodic boundary condition. For simplicity, we set the 
Planck constant h=l in this paper. 

First, we define a Hamiltonian with neighboring n- 
body interactions for n > 2 such as 

N 

h " : = £ ® • ■ ■ ® °f+n-l' ( 3 ) 

3=1 

where is a Pauli Z operator acting on the spin at the 

(n) 

j-th site, and Jj is a site-dependent coupling constant. 
Next, we define the most general Ising-type Hamiltonian 
H n for n > 2 composed of single-spin Hamiltonians and 
neighboring at most n-body interactions, 

n 

H n := H l + &mh m> (4) 

m=2 

where a constant A m denotes the strength of each h m 
for 2 < m < n and H\ := Y^jLi ^j a j denotes the contri- 
butions of the single-spin Hamiltonians. Since H± com- 
mutes with h m for all to > 2, it only generates additional 
local unitary operations and does not affect properties of 
entanglement. We set Hi = without loss of generality 
for investigating the entanglement evolution. 

Note that in physics, it is natural to expect that many- 
body interactions are weaker than few-body interactions. 
Thus, in many cases, we expect that a set of coefficients 
{A m } satisfies the relationship 1 > A m > A m > for 1 < 
to < to'. However, in special cases, this relationship is 
not necessary to be satisfied for all coefficients, since some 
of the coefficients have to vanish due to the symmetry of 
the system. For instance, for parity invariant systems, 
the odd terms h^k+i vanish. 

The state at time t is given by = e ~ lff "' |^ > 

where |^o) is an initial state. We study the entangle- 
ment distribution of the states ^(t)) by investigating 
the infinite-time average of the amount of entanglement, 
defined by 

1 f T 

(Emw) T;oo ~ lim - / E MW (\V{t)}), 

and the shortest time required for achieving (Emw)t-oo- 
We choose a separable initial state |^o) that generates 
(Emw)t-od as high as possible. In Appendix [Aj under 
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the assumption of phase crgodicity, where the phases of 
^(i)) are sufficiently randomized in the long-time limit, 
it is shown that 

( £ mw) T;co < 1 - 2lv^T' 

holds and the maximum value is achieved if and only if 
the initial state I'J'o) is given by ®£Li |+»), where 
are the eigenstates of with eigenvalue +1. Since this 
phase ergodicity assumption does not necessarily hold for 
H n in general, the maximum value is not guaranteed, but 
the initial state gives at least a lower bound for the gen- 
erated entanglement. Thus, we fix the initial state to be 
|^o) = ln [H| , it is shown that the maximum 

value of (Emw)t-oo ^ s g rea ter than the average amount 
of entanglement over the typical entanglement given by 
Eq. (TT]). That is, there is a possibility for the states 
to have an even larger amount of entanglement on aver- 
age than that of the typical entanglement. 

We denote the entanglement distribution of a set 
of randomly selected states {^(t))} generated by the 
Hamiltonian dynamics by {EmwW}(- We compare it 
with that of the typical entanglement of random states 
{^Mwjrand- We focus on the effects of the order of 
the many-body interactions, n, and the distribution of 
coupling constants {j[ m ^} for 2 < m < n. In order 
to examine their effects independently, we consider four 
Hamiltonians, h„, h n , H n and H n , where h n and H n 
are special cases of Hamiltonians of h n and H n with uni- 
form coupling constants. We denote the corresponding 
entanglement distributions by {E^ w (t)} t , {E^ w (t)} t , 
{E$b(t)} t and {E*fc(t)} t . 



III. THE CASE OF NEIGHBORING n-BODY 
INTERACTIONS 

In this section, we investigate the average and standard 
deviation of the entanglement distribution for the Hamil- 
tonian composed of neighboring n-body interactions h n 
defined by Eq. (J3J. We also show the time evolution of 
the two-point correlation functions for h n , which exhibit 
properties of correlations generated by the many-body 
interaction Hamiltonian. For simplicity, we drop the in- 
dex n of the coupling constant J- in Eq. ([3]) and denote 
the coupling constants by Jj in this section. 



A. Uniform coupling constants, h n 

We first investigate the case of the Hamiltonian h n , 
a special case of h n with uniform coupling constants, 
Ji = J. We denote the eigenstates of the Hamiltonian 
h n by |ai • • • cln) := ®£Li l a i): where the index a\ ■ ■ ■ ajy 
is a sequence of binary numbers a% G {0, 1} that cor- 
respond to the eigenvalues {— 1, +1} of erf. The basis 
{|oi ■ ■ • a at)} is usually referred to as the computational 
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FIG. 1. (Color online) £mw(|^(0)) generated by the uniform 
coupling constant Hamiltonian h n . The blue dotted-dashed 
line, black solid line and red dashed line indicate n — 2, 10, 20, 
respectively. It takes exactly t = 7r/4|J| to reach the maxi- 
mum, independently of n. 



basis. Then h n is written as 

N 

K = a f ® ■ • • ® a f+n-l 

3 = 1 

N j+n-1 

= J ^ X] ( 2a k- 1)|oi---ojv) (ai-"ffliv|- 

ai ,■ ■ ,«» 1=1 k=j 
=0,1 

From the initial state 

l*o)=^ 1 |+ l )=2- Ar / 2 ( 5 ) 

ai ,■ ■ ■ ,aw 

=0,1 

the time evolution of the state is described by 

|*(t))-2- w / 2 t[e- iJt nitr 1 ^-V\a 1 ...a N ). 

ai,— ,a N 3= 1 
=0,1 

Due to the translational invariance of the Hamiltonian 
h n and the initial state, all one-site reduced density ma- 
trices are identical. Thus, for any l-th spin, the reduced 
density matrix pi(t) = Tr^ ; |*(t)) (*(t)| is 

u . If 1 cos"2Jt\ 
pl{t > = 2 [cos n 2 Jt 1 )■ 

Then, we obtain 

£W(|*0)» = 1 - cos 2 ™ 2Jt. (6) 

Note that Emw (\^ (t))) is independent of the system size 
N. 

In this case, the entanglement evolution Emw (t))) 
for several different values of n is shown in Fig. [TJ It 
is seen that the maximum £mw = 1 is reached only at 
t = 7r/4|J| independently of n. However, larger values of 
n provide a faster speed to reach a neighborhood of the 
maximal value and a longer duration to stay in the the 
neighborhood of the maximal value. 
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As .EmwG^K^)}) i s periodic with period 7r/2|J|, it is 
easy to calculate the average amount of entanglement 



(Emw) t . 00 = 



1 



tt/2| J\ 



*/*\J\ jo 
1 (2n-l)H 
(2n)H 
1 



1 - cos 2 ™ 2Jtdt 



1 



where n!! = n(n — 2)(n — 4) ■ ■ • . The standard deviation 
ct-oo is obtained by 



0T;oo 



(4n-l)!! /(2n-l)!! 



(4n)! 



(2n)! 



1 1 



Because (EmwIt™ scales polynomially with n, the time 
evolution by h n cannot achieve the average amount of 
the typical entanglement even when n is large. 
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FIG. 2. (Color online) -E M w(|*(i)}) generated by site- 
dependent neighboring n-body interactions h n and for the 
number of spins iV = 50. The figures (a), (b), (c) and (d) 
correspond to n = 3, 4, 5, 6, respectively. The red constant 
line in each figure represents the upper bound of the infinite 
time average given by Eq. (|9]). The distribution of {Ji} is set 
to Gaussian centered at J with standard deviation 0.5J. 



B. Site-dependent coupling constants, h n 



Next, we investigate the Hamiltonian h n defined by 
Eq. in which the coupling constants J^ n ' arc site- 
dependent. Since the eigenstates of h n are the same as 
those of h n , the time evolution of the state from the initial 
state |$o) is given by 



A" 



air-- .aw j — 1 
=0,1 



The reduced density matrix pi (t) of the l-th spin is 



Pi(t) 



i ( i Mt) 



2 U(*i i r 



where 



A t (t)= cos2J fc t. 

fc=i-n+l 



(7) 



Note that p\ (t) now depends on I since the Hamiltonian 
h n is no longer translationally invariant. The entangle- 
ment evolution is 



1 N 

2WI*(*)>) = l-jy£(4(t)) 



(8) 



Although £'mw(|^'(0)) depends on the distribution of 
{Ji}, an upper bound of the long-time average is easily 
obtained by using 

1 f T n , fl (if = 0). 
Hm - / cosflidt - J 



T-*x> T 



(if 0^0). 



Then, the average amount of entanglement bounded by 

N 



(EMw) T ., 00 = l-^J2((Mt)) 2 ) 



>T;oo 



1=1 



1 -44E( II (l + cos4J fe t) 



where 



< 1 



A' 



2" AT 



2" ' 



i=l \fc=Z-n+l 



T;oo 



(9) 



IT (l + cos4J,t)\ >1 
;=1 \fc=i-n+l / T;oc 

is used to evaluate the last expression. For a completely 
random distribution of {Ji}, the equality holds. 

In Fig. [5J the entanglement evolution is shown when 
the distribution of {Ji} is Gaussian centered at J with 
standard deviation 0.5J. It is seen that entanglement 
concentrates around the value given by Eq. ([9]) . The long- 
time average and its standard deviation are numerically 
obtained by 



and 



(£mw) Too - 1 - 



0.07 



2 0.i 



(10) 



(11) 



The average amount of entanglement (Emw) t is greater 
than the upper bound given by Eq. @ for large n, but 
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this is merely an effect of finite time and it reaches the 
upper bound in the infinite time average. 

These values can be obtained in time of 0(1/ J) that 
is deduced from Eq. ((5]). By comparing Eq. (fTU|) with 
the average amount of the typical entanglement given by 
Eq. (JT]), they coincide if n ~ N. However, the standard 
deviation cr Too does not coincide with that of the typi- 
cal entanglement given by Eq. ©. We conclude that a 
set of states obtained by the Hamiltonian dynamics with 
h n can simulate the typical entanglement in terms of only 
the average amount of entanglement but not the standard 
deviation. Comparing to the average entanglement gen- 
erated by h n , this result indicates that the randomness 
introduced in the large n-body interactions is necessary 
to generate high average entanglement. 




FIG. 3. (Color online) The correlation function C x (r,t) as a 
function of the distance between the spins r and time t for the 
system evolving by the Hamiltonian h n and for N = 50 and 
n = 40. When t is small, two spins with large r is correlated. 
However, the correlation for large r rapidly decreases with 
time and disappears after a short time. 



C. Two-point correlation functions 



In this subsection, we investigate the dynamics of cor- 
relation functions, which are often discussed in associ- 
ation with entanglement in physics. In order to clarify 
their connection, we examine two-point spin correlation 
functions between the first spin and the (r + l)-th spin 
C w (r,t) for W £ X,Y,Z. Using a two-spin reduced 
density matrix of the first spin and the (r + l)-th spin 
denoted by pi_ r +i and a single spin density matrix of the 
i-th spin denoted by pi, C w (r,t) is written by 



C w (r,t) = Trp w aYa™ +1 - (Tr Pl af )(Trp r+1 a^ +1 ), 



(12) 

where the distance between the two spins r is only taken 
for r 
i ii t ion . 



for r = 1, • ■ ■ , y — 1 due to the periodic boundary con 



It is easy to see that for C x (r, t) and C Y (r, t), the first 
term of the right hand side of Eq. (fT2]l is a function of the 
elements (pi, r +l)i4, (pi, r +i)23, (pi,r+i)32 and (pi, r +l)4i 
in the a z basis. On the other hand, the Meyer- Wallach 
measure is a function of the diagonal elements (pi. r +i)u 
for i = 1, • • • ,4 and (pi,r+i)i3, (pi,r+i)24, (/>i,r+i)3i and 
(Pi.r+i)42- Thus, C x (r,t) and C Y (r, t) reflect aspects 
of correlations of the system not included in the Meyer- 
Wallach measure. On the other hand, Trpi !r+ icrf cr z +1 is 
a function of the diagonal elements, so it is related to the 
Meyer- Wallach measure. 



For the system evolving by the Hamiltonian h n , the 
two-point spin correlation functions are analytically cal- 
culated to 



C z (r,t) = 0, 
C Y (r,t) = 



and 
C X (r,t) 

^ cos 2 Jit — cos 2 Jit cos 2 Jyt (if r < n), 







leP(N) 

l'eP(r) 



(if r > n), 
(13) 



where 



L = 



' [N - n + 2, N - n + r + 1] U [2, r + 1] 

(if r < N-n + 1), 
[N - n + r + 2, N - n + r + 1] U [2, N - n + 1} 

(if r >N-n+l). 

and P(k) = [k — n + 2, k + 1]. Except for the X direc- 
tion, distant spins cannot be correlated by the Ising-type 
Hamiltonian h n . 

In Fig. [3] and Fig. |4l we show how C x (r,t) decreases 
with r and t, depending on the order of interaction n. 
When n is small, C x (r, t) oscillates for small r and is zero 
for large r, which is implied by Eq. (fT3"|) . On the other 
hand, as shown in Fig. [3l when n is large and t is small, 
C x (r, t) is non-zero even for large r . Thus the large 
neighboring n-body interaction can create strong corre- 
lations between distant spins, as well as a large amount 
of entanglement. However, in contrast to the entangle- 
ment evolution, the correlation rapidly decreases with 
time and the correlation between distant spins vanishes 
after a short time. Thus, applying the large neighbor- 
ing n-body interaction does not dramatically enhance the 
two-point spin correlation functions in this model. 



IV. THE CASE OF UP TO NEIGHBORING 
n-BODY INTERACTIONS 

In this section, we investigate the average and standard 
deviation of the entanglement distribution for the general 
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C A (r,t) 
1.0 




FIG. 4. (Color online) Correlation functions C x (r) as a func- 
tion of the distance r for a system of N = 50 spins. The 
figures (a), (b), (c) and (d) show n = 45,30,20,5, respec- 
tively. In each figure, the diamond, square and circle points 
represent t = 7r/(50J), 7r/(30J), 7r/(10J), respectively. It is 
seen that, for r > n, the correlation function is strictly zero 
in all cases, and that the correlation decreases quickly with 
time. 



Hamiltonian H n , composed of up to neighboring n-body 
interactions h n , denned by Eq. (j4}. 



A. Uniform coupling constants, H n 

We first investigate the case of the Hamiltonian H n , a 
special case of H n composed of up to neighboring n-body 
interactions with uniform coupling constants Jj = J- 
The Hamiltonian is given by 

n 

H n = A m h m . (14) 

m=2 

Since the Hamiltonian H n is translationally invariant for 
uniform coupling constants, the reduced density matrix 
pi(t) does not depend on I and is in the form of 



P(t) 



1 / 1 B(t) 



2 \B(ty 1 



where B(t) is given by 

n 

B{t)=2^ y, n 



Oiv ,ojv m=2 
=0,1 

n j+n— 1 

x Yi exp[-2zA m Y[ (2o fc -l)Jt], 

j=n—m+l k=j 

and B* is the complex conjugate of B. The entanglement 
evolution is 

£mw(|*W)) = i-|bWI 2 , 
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FIG. 5. (Color online) £mw(|*(£)}) generated by H„. The 
figures (a), (b), (c) and (d) correspond to n = 3, 4, 5, 6, respec- 
tively. The {A m } are set to 2e/m with e = y/3/10 and the 
figures are plotted up to t = (27rn 4 )/e which is 0(poly(A^ 1 )). 



which is independent of the number of spins N. 

Although the average amount of entanglement is dif- 
ficult to calculate in this case, upper and lower bounds 
of the average amount of entanglement are analytically 
obtained by 



1 - 2 "("-l) " ^ MW ^;oo 



< 1 



2 2(«-l)' 



(15) 



where a = log 2 | ~ 1.42 We show the derivation of these 
bounds in Appendix. [EJ Similarly, the standard devia- 
tion is estimated by 



0T;< 



0{2- 2n ) 



These values do not depend on the details of the strengths 
each h m and A m in H n . In Fig. [5l the entanglement evo- 
lution generated by H n with A m = l/poly(m) is plotted 
for Jt £ [0, 0(l/poly(A„))], along with its time-average. 
The entanglement evolution generated by the Hamilto- 
nian with exponentially decreasing A m behaves similarly 
to the case with polynomially decreasing A m . 

With numerical calculations, we obtain the average 
and standard deviation of the entanglement distribution 



2&(n-l) ' 



and 
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2<5(n-l) 



where 



(a J) 



(2.0,1.4) (if A, 
(1.9,1.6) (if A, 



l/poly(m)), 
1/ exp(m)), 



(16) 



(17) 



(18) 
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and 



i\lW 
] _ 



Hi) 



(1.8.1.4) (if A, 

(1.7.5.5) (if A, 



l/poly(m)), 
1/ exp(m)). 



These averages are taken over < t < T n , where T n is 
chosen to be a polynomial of (JA ra ) _1 . Eq. (fT8| shows 
that the average amount of entanglement is not strongly 
dependent on {A m }, which is expected from the bounds 
given by Eq. (fT5"j) . The average (£mw) t can be as high 
as that of the typical entanglement if n ~ N/a, namely, 
n ~ N/2. In that case, the standard deviations are com- 
parable with that of the typical entanglement as well. 
Hence, if n is set to ~ N/2, the entanglement distri- 
bution obtained by the time evolution generated by H n 
simulates that of the typical entanglement in terms of 
both the average and the standard deviation. 

In order to investigate the time Too required to reach 
(-E'Mw)t-oo' we cons ider the time evolution operator 
U n = e~ tHnt which can be decomposed into 



e -iA 2 h 2 t e - 



,— iA n hnt 



It is obvious that 



,-jA./i„( 



is almost the identity, if 



A n t is negligibly small compared to the largest eigen- 
value of h n , which is equal to NJ. Hence, it takes at 
least, 0(l/(NJA n )) to achieve (Smw)too- The nu- 
merical results shown in Fig. [5] support the claim that 
0(l/(Jpoly(A„))) is enough to reach the infinite-time 
average. For n ~ N/2, if A m scales polynomially with 
m, the time is a polynomial in N. 



B. Site-dependent coupling constants, H n 

We show the result of the most general Hamiltonian 
H n defined by Eq. In this case, each h m contained 
in H n has coupling constants depending on the site i, 

J- m) - £mw(| *(*))) is obtained by 

iV 



^mw(i*w)) = i-Ei^wi 2 ' 

(=1 

where Bi(t) is given by 

n 

Bl (t)=2^ yi n 

ai,— ,a N m =2 
=0,1 

I j+n-l 

x J] exp[-2iA m _ 2 Y[ (2a fc -l)j/ m) i]. 

j—l — ra+l k=j 

The entanglement evolution as well as the average is 
shown in Fig. |5J Numerically we obtain 



(%w) r — 1 



2 2.12(n-l) ! 



and 
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FIG. 6. (Color online) £mw(| *(£))) generated by H„. The 
figures (a), (b), (c) and (d) correspond to n = 3,4,5,6, 
respectively. The plot is shown for the time scale t € 
[0, poly(l/(JA n ))]. The distribution of \Ji} is given by a 
Gaussian centered at J with standard deviation 0.5J. The 
{A m } are set to 2e/m with e = y^/lO and the figures are 
plotted up to t = (27rn 4 )/e which is 0(poly(A^ 1 )). Compar- 
ing to Fig. \5\ the average amount of entanglement is higher. 



The strengths of the many-body interactions are cho- 
sen to be A rn = 2e/m where e = y/3/10 and the inte- 
gral interval is taken from to Jt = 27m 4 /e which is 
0(l/(Jpoly(A n ))). 

As is expected from the results of the previous sub- 
sections, both the average and the standard deviation 
of the entanglement distribution coincide to those of the 
typical entanglement by choosing n ~ N/2. The time 
required for achieving a neighborhood of maximal en- 
tanglement is deduced from the same argument in Sub- 
section IIV Al and numerical calculations confirm that 
= 0(l/(Jpoly(A„))) is enough to achieve the infinite- 
time average if {A„} scales polynomially with n. That is, 
(Emw)t-oo ana - a T;<x> are obtained in a time polynomial 
in N. 



V. ENTANGLEMENT DISTRIBUTION 

In this section, we compare the entanglement distribu- 
tions obtained by Hamiltonians h n , h n , H n and H n and 
that of the typical entanglement of random states. In 
Table. HI we summarize the results obtained in the previ- 
ous two sections. First, we note that for all Hamiltonians 
we analyzed, the obtained entanglement distributions de- 
pend on n but not on N. 

For the average of the entanglement distributions ob- 
tained by the Hamiltonians, it is shown that the Hamilto- 
nians h n , H n and H n can generate entanglement as high 
as that of the typical entanglement by choosing suitable 
n. For h n , n ^ N is required for the average to be com- 
parable with (j5Mw) rall( j- For H n and H n , n ~ N/2 is 
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TABLE I. The time-averaged entanglement for each Hamilto- 
nian defined in Section fll Bl {Emw) t . co is an approximation, 
found analytically, of the average in the infinite time limit. 
Numerical results are given by (E M w) T and a Too , which de- 
note an average over the time interval to Too and its stan- 
dard deviation, respectively. The time required to reach the 
infinite-time average (Emw} t . <x is denoted by Too. 



(£mw) 5 



(£mw) To 



&t 



It ,1 



1 - 



1.7 



0.07 

20. 6n 



1 



2«(n-l) 



(a€ [log 2 |,2]) 



1.4 



22.0(n-l) 



1.4 



21.8(n-l) 



1.1 



22.1(7!-!) 



1.1 



21.7(n-l) 



sufficient for achieving (Emw),^. On the other hand, 
for the Hamiltonian composed of constant neighboring 
n-body interactions h n , the infinite-time average entan- 
glement (-Emw)t'oo sca l es polynomially with n even when 
n is set to its maximum value N. 

For the standard deviation of the entanglement distri- 
bution, the results for h n and h n show that neighboring 
n-body interactions alone do not simulate the standard 
deviation of the typical entanglement. For H n and H n , 
n ~ N/2 is also sufficient for achieving the standard de- 
viation coinciding to that of the typical entanglement of 
random states <7 r and- Thus these Hamiltonian dynamics 
can provide a set of states simulating the entanglement 
distribution of {-EMw}rand hi terms of up to the second 
order of distribution, namely, the average and standard 
deviation. 

However, these results do not guarantee that the entan- 
glement distributions themselves obtained by the Hamil- 
tonians coincide with that of the typical entanglement. 
We show the calculations of the entanglement distribu- 
tions for N = 8 for h n , H n and the typical entanglement 
in Fig. [7] For each Hamiltonian, n is set to the value 
that achieves the average amount of the typical entangle- 
ment. The entanglement distributions obtained by these 
Hamiltonians are clearly different from that of the typi- 
cal entanglement for N = 8. In these cases, the entangle- 
ment distributions obtained by Hamiltonian dynamics is 
concentrated around the values higher than that of the 
distribution of the typical entanglement. Such properties 
of the entanglement distributions obtained by Hamilto- 
nian dynamics may be useful for applications requiring 
to use a set of states with a higher concentration of en- 
tanglement in terms of Meyer- Wallach measure. 

Next, we analyze the shortest time t^ required to 
achieve (Bmw)t-oo' ^or ^« ano - h n , T oc is 0(1/ J). For 
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FIG. 7. (Color online) The entanglement distributions for 
N — 8. The sampling number is 10000. The figure (a) is the 
distribution obtained from the typical entanglement of ran- 
dom states, {i?Mw}rand- The figure (b) is the distribution 
obtained by the Hamiltonian dynamics of h„, {-Ejyjyy- (*)}t, 
where n = 6. The figure (c) is the distribution obtained by 
the Hamiltonian dynamics of H„, (B^fi)}!, where n = 4. 
Although the lower-order moments coincide with those of the 
typical entanglement, as shown in Table. [I] the entanglement 
distributions are different, which implies that the higher-order 
moments must be different. 



h ni J is the center of Gaussian of the distribution of 
the coupling constant {Ji}- Thus the average amount 
of entanglement is achievable in a time scales indepen- 
dent of N. On the other hand, for H n and H n , is 
inversely proportional to the strength of the neighboring 
n-body interaction, A„, which is determined by the de- 
tails of the Hamiltonian. If A m scales with l/poly(m), 
then Too = poly(n). Because n ~ N/2 is necessary in 
order to achieve (Emw) rand' ^ ne time required to achieve 
(-Emw)t oo i s a ^ so polynomial in N. On the other hand, 
if A m scales with l/exp(m), the time required to achieve 
(Emw)t-oo by the Hamiltonian dynamics is exponential 
in N. 



VI. SUMMARY 

In this paper, we have studied the entanglement dis- 
tributions obtained by one dimensional spin- 1/2 Ising- 
type Hamiltonians composed of several types of many- 
body interactions. We have shown that, when the time- 
independent Hamiltonian is composed of up to n ~ N/2 
neighboring ?i-body interactions, a set of states ran- 
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domly selected during the time evolution generated by 
the Hamiltonian can simulate the entanglement distri- 
bution of the typical entanglement of random states in 
terms of both the average and the standard deviation. 
Our results imply that n ~ logiV neighboring n-body 
Ising-type interactions are not sufficient for generating 
high entanglement, even though the order of many-body 
interactions n scales with the number of spins N and the 
interactions are not local interactions in the sense defined 
by [IH]. We have also shown that the time required to 
achieve such a distribution is polynomial in the systems 
size if the strength of neighboring m-body interactions 
scales with l/poly(m). 

On the other hand, when the Hamiltonian is composed 
of only neighboring n-body interactions, site-dependent 
coupling constants and n ~ N neighboring 71-body inter- 
actions are required to simulate the average amount of 
entanglement of the typical entanglement. Although the 
standard deviation does not coincide to that of the distri- 
bution of the typical entanglement, the time required to 
achieve the distribution to simulate the average amount 
of entanglement alone is shown to be independent of the 
system size in this case. 
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Appendix A: Derivation of the initial state 

In this appendix, we show that the initial state \^o) = 
®£Li l+i) gives the maximal average amount of entangle- 
ment (Emw)t -oo' ^ the eigenstates of the Hamiltonian 
and the initial state are all separable. 

In (2ll | , we have introduced phase-random states, a set 

of pure states \%2 a =i r a& %Va \ u a)}{p a } w hh fixed ampli- 
tudes r a > and uniformly distributed phases {<p a } in 
a fixed basis {|m q )} and derive a formula for the average 
amount of linear entropy of entanglement of phase ran- 
dom states (i?L) p haso for every possible bipartite parti- 
tion. We have analyzed the entanglement evolution gen- 
erated by a general Hamiltonian denoted by 

H = ^2e a \e a ) (e a \ , 

a=l 

where {|e Q )} are eigenstates, from an initial state denoted 



by 

\*{0)) = Y,r a e**> \e a ) , (Al) 

a=l 

where r a ,uj a £ R and |r a | 2 = 1. Under the assump- 
tion of phase ergodicity in the sense that the distribution 
of phases of e"' £ ° t+ ™° are uniform in [0, 2ii\ in the long- 
time limit, the infinite-time average of the amount of 
entanglement (El)t-,oo coincides to the average amount 
of entanglement of phase random states (i?L) p haso- From 
the formula for (El) phase, it is straightforward to derive 
the following lemma in terms of the Mayer- Wallach mea- 
sure of entanglement. 

Lemma 1 If the system exhibits phase ergodicity, then 

(Emw)t:oo = (^MW) p hasc 

N 2 N 

= ^E[^(/ 3 av) + ^(Pa V fc )-^(Pav)]-E r ^Mw(|£ a )). 
k=l o=l 

(A2) 

where 

2 n 

a=l 

pl v := Tr-.fepav and := Tr fcj o av for k =!,-■■ ,N. 

We use this lemma for the Ising-type Hamiltonians and 
for separable initial states. We denote the eigenstates of 
the Ising-type Hamiltonian H n in the Pauli Z basis by 
using the binary representation \a) defined by 

I a) := |ai • • • a N ) . (A3) 

Using this notation, the time evolution of the state |\&(t)) 
generated by H n from a general initial state given by 
Eq. (|A1|) is written by 

2 N 

m))=J2r a e i ^- E ^\a), (A4) 

a=l 

where E a is the eigenenergy corresponding to \a). 

Assuming phase ergodicity, the average amount of en- 
tanglement (Emw)t-oo * s obtained by using the lemma. 
Since the eigenstate \a) is a separable state for all a, 
EMw(\a)) = 0, (£mw) T;(X) is given by 

2 N 

(Emw)t:oo = m Zj^Pav) + ^(Pav ) - St(Pav)]. 
k=l 

(A5) 

where 

2" 

Pav = l«) ("I ' 

a=l 
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The right hand side of Eq. (IA5|) can be written in terms 
of r a by 



JV 



(Emw) t . i00 = ^£[£ 



r Pk(i,m) ' Sfc(j,n) 



k— 1 (i,m) 



EJ r 2 1 

p k (i,m) s k (i,m)l 

(i,m) 



(A6) 



where pk(i,m) := i + 2m ■ 2 N k and Sk(i,m) := i + 
(2m + 1) ■ 2 N ~ k and the summations are taken over i,j = 
I,-- - ,2 N - k and m, n = 0,---2 k ~ 1 - 1. By using the 
normalization condition of r a , the first term of Eq. (|A6[) 
is evaluated as 

— V V r 2 r 2 
jyZj Pk{i,m)' s k (j,n) 

k—1 (i,m) 

Ci.n) 

4 W 



Pk(j,n) 



fe=l (i,m) 



< 1. 



1/2 



The equality holds if and only if r 2 fc(jm) 
for any 

On the other hand, a lower bound of the second terms 
of Eq. (|A6[) is given by 



— "V "V r 2 r 2 > 4 

/ i / t pk(i,m)' Sfe(i,m) — ^ 

/c— 1 (i,m) 



JV 



JV 



fc— 1 (i,m) 



r pk(i,m)'sk(i,m.) 



since the arithmetic mean is always greater than or equal 
to the geometric mean. The equality holds if and only 

if R(r) := Y,{i, m ) r l k {i, m ) r2 s k {i, m ) does not de P cnd on k - 
By substituting these relations, the upper bound of the 
average is obtained by 

(Emw) T:oc < 1 - 4i?(r), 
where the equality holds if and only if the two conditions, 



Pk{i,nn) 



1/2, 



(i,m) 



and 



(i,m) 



fc(z,m)' Sk(i,rn) 



R(t), 



(A7) 



(A8) 



are simultaneously satisfied. 

Since we consider a separable initial state 

= Ea=i r a e!IJ ° l a ), ^ must satisfy the con- 
dition rankTr^ fc |*(0)) (*(0)| = 1 for any k, that 
is, 

Vfc.rank^ - ^ = 1, (A9) 



where Cl k (i,m) := r Pk ^ !m )r Sk ^ m )e t( - w "^^ w «k«.'»>) and 
m) is the complex conjugate of £lk{h m )- Under the 
conditions given by Eq. (|A7|) and Eq. (|A9| , we minimize 
i?(r) so that we can obtain the maximum value of the 
average (Emw)t-oo f° r separable initial states. 
From Eq. (|A"7|) and Eq. (JX9J, we obtain 

| 5^fifc(i,m)| = i 

With this condition, i?(r) = ^(i,m) l^fc(*, m )| 2 i s easily 
maximized by Lagrange's method of undetermined mul- 
tipliers such as 



R(r)> 



1 



2JV+1 ■ 



where the equality holds if and only if |fifc(i,m)| = 1/2 N 
and Argr^fc(i,m) is constant for any k, i and to. There- 
fore, we have proven that, for any separable initial states, 



(£mw) T ;oo ^ 1 



1 



)JV-1 1 



The equality holds if and only if £(i,m) r l k {i, m ) = 1 / 2 
for any k, and |f2fc(i, m)| = 1/2 W for any k, i and to, 
which are equivalent to r a = 1/2 N ' 2 and lu ci = lo for any 
a. Thus, the initial state |* ) = e iuj <g)f =1 gives the 
maximum amount of average of entanglement. By drop- 
ping the global phase e lw , we obtain |*o) = ®£Li !+<)• 



Appendix B: Upper and lower bounds of the average 
entanglement for H n 

In this appendix, we show the following statement; for 
the Hamiltonian H n defined by Eq. 

4 4 

1 " oa(n-l) " ( £; Mw) T ;oo < 1 " 9 2(„-l)> 



where a = log 2 | ~ 1.42. 

From the initial state |^o) given by Eq. ([5]), the time 
evolution of the state by the Hamiltonian H n is formally 
written by 



(Bl) 



where \a) is the binary representation of the eigenstates 
of H n (and also H n ) introduced by Eq. (|A3j) . and E a is 
the cigenencrgy corresponding to \a) given by 

n JV j+m— 1 

£ a = J^A ro ^ J] (2a t -l). (B2) 

m=2 j'=l t=j 

Since the Hamiltonian -ff„ is translationally invari- 
ant, the Meyer- Wallach measure £'mw(|^ , (0)) i s given 
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by E MW (\^(t))) = 2 - 2Trp 2 (i). By calculating the re- 
duced density matrix pi(t) of \$(t)) given by Eq. (|B1I) . 
the Mayer- Wallach measure is obtained by 

3mw(|*(*)» = 1 - C0S [( £ a ~ £&)*], (B3) 



,6=1 



where e a = E a — E a+2 N-i. From Eq. (|B2|) . e a is given by 

n AT+1 JV j+m-1 

m=2 j=N-m+2 t=j s=2 

(B4) 

To calculate the bounds of (£mw)t-oo' ^ ^ s sufficient 
to evaluate the term (cos[(e/t — ^i)t]) T . 00 . This term gives 
one for e/. = ei and zero for 7^ e;. By defining a 
quantity S by 

S := (the number of (fc, /) such that = ei ), 

we obtain 



(Emw) T;00 - 1 - 2 2(JV-1)' 



From Eq. (|B4[) . we see that e a depends on a t for 
t = N-n + 2,N-n + 3,---,N-n + n where 
n = 1,--'JV, but it does not depend on a t for t = 



1, n + 1, ■ ■ ■ , iV — 71+ 1. By introducing a vector notation 

Q(iV-Tt+2 )-m : ~ { a N-n+2,' ' ' , a N, a 2i''' ,«ti), e a — e 6 

in Eq. (|B3j) is only determined by a(jv— n+2)-m ancl 
&( AT-n+2)->n- We define a quantity £, which is the number 
of pairs (a (N _ n+2) ^ n , b {N _ n+2) ^ n ) satisfying e a = e b . 
Then c can be expressed in terms of ^ by 



5 = £ x 2 2 ( Ar - 2n + 2 ), 

where the factor 2 2 ( 7V_2n + 2 ) appears due to the choice of 
at and b t for t = 1, n+1, • • • , 2V— Recall that e Q (e^) 
is independent of a t (bt) for t = l,n + 1, • ■ ■ ,N — n+1. 
Thus, we obtain 



(E-. 



MW, 



T;oo 



1 



24n— 6 ' 



We evaluate the upper and lower bounds of £. For the 
lower bound, when a = b, e a = is trivially satisfied. 
Since the number of choices of a(jv_ n +2)-m is 2 2 ^"~ 1 \ we 
obtain 2 2 (™~ 1 ' < £. It is cumbersome to derive the upper 
bound of £, but direct investigation of the expression of 
Eq. JB4| provides a bound £ < 6 n_1 . Therefore, the 
upper and lower bounds of (-Emw)t-oo are 



n-l 



< (Emw) 



T;oo 



< 1 -4 



n-l 
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